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Abstract 

It has been proven by Rosu and Cornejo-Perez [1] [2] that for some nonlinear second-order ODEs it is a 
very simple task to find one particular solution once the nonlinear equation is factorized with the use of 
two first-order differential operators. Here, it is shown that an interesting class of parametric solutions 
are easy to obtain if the proposed factorization has a particular form, which happily turns out to be 
the case in many problems of physical interest. The method that we exemplify with a few explicitly 
solved cases consists in using the general solution of the Riccati equation, which contributes with one 
parameter to this class of parametric solutions. For these nonlinear cases, the Riccati parameter serves 
as a 'growth' parameter from the trivial null solution up to the particular solution found through the 
factorization procedure. 

PACS: 02.30.Jr, 02.30.Hq, 11.30.Pb 



1 Introduction 

Despite powerful integrability methods, such as the Lie group-theoretical approach, Painleve analysis, exis- 
tence of Lax representations and the associated inverse scattering transforms, the task of obtaining solutions 
of nonlinear second order partial and ordinary differential equations (ODEs) remains one of the most difficult 
problems in mathematical physics; in some cases, even finding one particular solution turns out to be a very 
difficult matter [3113]. However, in a number of cases, it has been proven that finding one particular solution 
turns out to be easier than expected. In the case of polynomial non-linearities, Rosu and Cornejo-Perez 
[1] [2] , working with a factorization procedure stemming from a work of Berkovich [|5i] , have found that if the 
second order nonlinear differential equation can be factorized into two first order differential operators then 
it is easy to find a first particular solution for the problem. They considered a nonlinear equation of the type 

u + g{u)u + F{u) ^0 , (1) 

where the dots represent derivatives with respect to the independent variable r, which is usually the traveling 
coordinate of a reaction-diffiision equation [T| . The method they proposed was to factorize this equation in 
the following form 

[D-Mu)][D^Mu)]u^O , (2) 
(where D = ^) which implies the following conditions on the functions (l)i{u) 

' ^+-^2 + ^^^) - gin) (3) 

<^102 = (4) 
u 

If Eq. Il]) can be factorized as in Eq. ([2]), then a first particular solution, say ui, can be easily found by 
solving 

[Z?-(/.i(m)]u = . (5) 
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2 Riccati- parameter solutions 



Of course, obtaining one solution of a nonlinear second order ODE does not guarantee that one may find 
more general solutions. However, what Rosu and Cornejo-Perez found was that in many cases (some of which 
will be described below) the function 0i(it) turned out to be a linear function of the dependent variable u. 
Hence, Eq. ([5]) turns out to be a Riccati equation for this variable, which is very fortunate, since we already 
know how to find the general solution for this equation once a particular solution is known. 

The appearance of the Riccati equation in linear second order differential equations is very common. In 
particular, it was very successfully exploited by Mielnik to find potentials which are isospectral to the simple 
harmonic oscillator potential [6], and it is a cornerstone for all SUSY developments [7]. However, it has not 
been used to solve nonlinear second-order differential equations at least in the way we present here. 

Thus, if 01 is of the form 4'i{u) = ciu + C2, Eq. ^ transforms into the Riccati equation 

ii — ciu^ — C2U = (6) 
and if a particular solution ui of this equation is known, then the general solution can be found as [B] 

ux,ci = "1 + T r , (7) 

where ^ 

/l(r) EE / {2ciUl{T')+C2)dT' (8) 



To 



and 

/2(t)= r e'^'-^'^dr' . (9) 

J To 

Notice that for 

Xs = cihir) , (10) 

a singularity may develop. 

Eq. ([7]) provides in turn what we call as a Riccati-parameter solution of the nonlinear equation. We 
notice that the first parameter, ci, is essentially the slope of the factorization function 0i, whereas the A 
parameter can be chosen in such a way as to prevent this solution from possessing singularities [B] , although 
for nonlinear differential equations this is not an absolutely prohibitive issue. It will be seen in the examples 
given in the following that the latter parameter acts like a label in this class of solutions placing them 
between the trivial null solution and the particular solution given by Eq. 



3 Examples of physical interest 

In this section wc find the explicit form of the Riccati-parameter solution for three nonlinear equations of 
physical interest that are polynomial type Lienard equations, i.e., similar to Eq. ([T]) but with F{u) a poly- 
nomial of order two and three in our cases. 

3.1 Modified Emden equation 

We start with the modified Emden equation 

■it + a u M + /3 u'"^ = , (11) 

for which the first rigorous study has been done by Painleve more than a century ago [8] who got solutions 
for P — /9 and f3 — ~o? . Recently, Chandrasekar et al |9j provided a detailed discussion of this equation 
from the point of view of the modified Prelle-Singer procedure that gives the construction of the solution 
in terms of elementary functions if such a solution exists [10], although lacono [11 noticed that much 
simpler connections with the Abel equation could be used to get the solutions. For the remarkable physical 
applications, see [9]. 

Employing (f)i = aiy/Pu and 02 — ai^y/j3u, where ai — — "^^^ , one particular solution is [2] 

"1 = tJ" r • (12) 

aiVP(T - To) 



2 



Hence, by using Eq. ([7]) one can find that the two-parameter solution of Eq. (fTTj) is 

1 1 



aiV^(T-To) A(7 



(13) 



In this case it is instructive to notice that when |A| runs from zero to infinity, the u\ solution goes from the 
trivial solution u = to the particular solution u = ui, as can be deduced from Eq. (|13|) and graphically 
seen in figure [TJ 




Fig. 1: ux{t) in the case of the modified Emden equation for A — —0.2, —0.4, —1, —6, from top to bottom respectively, 
and ai ~ —1. 

3.2 Convective Fisher equation 

We pass now to the convective Fisher equation [12] 

il + 2{iy - ^u)u + 2u{l - u) = . (14) 

The second term corresponding to convection is introduced to describe mechanical transport in competition 
with diffusive transport or cases when external bias fields are present. In the context of population dynamics 
which is typical for the Fisher equation, it has been used by Walsh et al [13j to simulate the population 
mobility according to spatial gradients in the food supply. 

Rosu and Cornejo-Perez found that for v — fj,/2 + the factorization functions are (pi = — /i(l — u) 
and (j)2 — ~2/^. Hence, a particular solution for this equation will be [2] 



Ml = [1 ± exp (/^(r - To))] ^ 
The A-parameter solution can be readily obtained in the form 



|-g-M(r~ro) ± 1] (^g-t^{T-To) ± l) - l] 

Once again, uq = and Uoo = ui. For the graphical representation see figure [21 



(15) 



(16) 



3.3 Generalized Lienard equation 

Consider now the generalized Lienard equation for cubic nonlinear oscillators 

zi + g[u) u + F3 , 



(17) 
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Fig. 2: Convective Fisher u\{t) for A = —0.2, —0.4, —1, —6, from top to bottom respectively, and fj, = 1. 
where F^^u) = Au + Bv? + Cv? . The previous equations can be seen as particular cases of this one. With 



A = \/B'^ — AAC, Rosu and Cornejo-Perez found that using 

'B + A 



Ci 



-i(B-A 



for 



B + A B-A _i ,^ 



(18) 



the following particular solution could be obtained 2J 
ui(t) 



B+A 
2 



B+A 



cxp(ai^(T-To)) 



exp (-ai^(T - To)) - C 1 - Cexp (ai^(T - To)) ' 
Now, using Eq. ([7]) and denoting ^ = r — tq, we can see that the two-parameter solution in this case is 

^exp(-ai^C) 



uxiO =ui + 



[exp (-ai^C) - C] [(A^ - l) exp (-ai^^) - AC^] 



(19) 



(20) 



Varying A between and oo, the Riccati-parameter solutions go from the null solution to ui. Plots of u\ for 
several values of A are displayed in figure [31 



3.^ Other cases of physical interest 

The examples we have provided here are not the only possible cases that can be solved with this method, 
but they show the typical solutions to be found. Other examples where Riccati-parameter solutions can 
be obtained in this way are: the generalized Burgers-Huxley equation with 5 1, 7 = 1, the generalized 
Fisher equation, with n = 2, the Dixon- Tuszynski-Otwinowski type equation, with n = A, and the Fitzhugh- 
Nagumo equation [T]. The linear factorization functions (pi of all these cases are given in [Tj. Last but not 
least, we would like to comment on the possible physical interpretation of the Riccati parameter A. We 
follow the works of Barton et al |14j and Monthus et al to assert that A is related to the introduction of 
finite interval boundaries on the abscissa of the problem. This is very well described in Section II A of [15] 
to which the interested reader is directed. Essentially, the introduction of boundary conditions at certain 
points on the axis generates a modulation of the particular solution as presented here and the A parameter 
can be fixed through the boundary conditions. When the boundary is sent to infinity the original particular 
solution is recovered. 
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Fig. 3: Cubic Lienard u\{(^) solutions for A = 0.7, 0.9, 1.1, 1.3, from top to bottom respectively, in the case A — C — 1, 
-B = 2 (A = 0) and ai = 1. 

In conclusion, we introduced here an interesting class of parametric solutions of a number of physically 
relevant nonlinear differential equations. They cover the space between the null or constant solution and the 
particular solution obtained by a simple factorization method proposed previously. 

We would like to thank Dr. O. Cornejo-Perez for a careful reading of the first draft of this work. The 
second author wishes to thank CONACyT for partial support through project 46980. 
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